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Abstract. Consider two manifolds A/"" and A'^" and a first-order Lagrangian L{u) 
for mappings u : M ^ N, i.e., L is an expression involving u and its first derivatives 
whose value is an m-form (or more generally, an m-density) on M. One is usually 
interested in describing the extrema of the functional C{u) = Jj^jL{u), and these 
are characterized locally as the solutions of the Euler-Lagrange equation Ej^ (u) = 
associated to L. In this note I will discuss three problems which can be understood as 
trying to determine how many solutions exist to the Euler-Lagrange equation which 
also satisfy L{u) = <&, where <E> is a specified m-form or m-density on M. The first 
problem, which is solved completely, is to determine when two minimal graphs over 
a domain in the plane can induce the same area form without merely differing by a 
vertical translation or reflection. The second problem, described more fully below, 
arose in Professor Calabi's study of extremal isosystolic metrics on surfaces. The 
third problem, also solved completely, is to determine the (local) harmonic maps 
between spheres which have constant energy density. 



0. Introduction 

In the study of pde, the search for solutions with special properties has a long 
history. Adding extra conditions, such as symmetry assumptions, frequently re- 
duces the complexity of a problem to a point where it can be explicitly solved. The 
resulting special solutions are then often used as the basis for formulating conjec- 
tures and as comparisons to the general solution so as to gain information about 
global behavior, asymptotics, long time existence, etc. In another direction, one 
is often interested in uniqueness questions, i.e., determining how much extra infor- 
mation must be specified in order to determine a unique solution to the pde under 
study. 
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In this note a combination of these two approaches is considered. I am going to 
describe three problems in which the PDE under consideration is the Euler-Lagrange 
equation for some first order functional and the extra condition takes the form of 
specifying the value of the Lagrangian density of the solution in advance. Normally 
this extra condition so overdetermines the problem that there will be no solutions 
at all. What I will be interested in will be the special cases for which, not only do 
solutions exist, but uniqueness fails. 

Let me describe a simple, concrete case of the general problem. Suppose given 
a first order functional of the form 



(1) 



jC{z) = / L{x,y,z,Zx,Zy)dx hdy 

J-D 



where "D C is a plane domain and z is & (real-valued) function on this domain. 
For any given function F on V, the problem is to describe the space of solutions to 
the Euler-Lagrange equation for C which also satisfy the first order condition 

(2) L{x,y,z,z^,Zy) = F{x,y). 

Some basic insight into the nature of this problem can be seen as follows. Differ- 
entiating (2) with respect to x and y yields the equations 

~t~ ~l~ ^xx "I" ^xy ~ ^x 

Ly -\- Zy ~\~ Lp Zxy ~\~ Lq Zyy Fy 

while the Euler-Lagrange equation of (1) is of the form 

(4) Lpp Zxx ~t~ '^Lpq Zxy ~t~ Lqq Zyy — L 

where L only involves x, y, z, p = z^ and q = Zy. The three equations (3-4) can 
be solved for the second partials (i.e., the 2-jet) of z in terms of the lower order 
partials (i.e., the 1-jet) of z provided that the obvious determinant does not vanish 
(which will, in general, be the case). The usual methods of total differential equa- 
tions [BCG] then imply that there is at most a three-parameter family of solutions 
to these three equations. Coupled with the original first-order equation (2), this 
results in a maximum of two parameters worth of solutions to the system. 

However, for the generic choice of F, the system (2-4) will be incompatible and 
will have no solutions at all. The question is to determine those F for which this 
system is sufficiently compatible that there will exist more than one solution. The 
"extra" compatibility of this system will turn out to be expressed as a system of 
partial differential equations on _F, which need not have any solutions. 

The first two problems to be treated below were suggested to me by Professor 
Calabi and are precisely of the type described above. In the first, the Lagrangian 
is the area Lagrangian for a graph z = u{x, y) over a plane domain, i.e., 



L = s/l 



Thus, the problem is to determine when two minimal graphs can determine the same 
area function. Obviously, two minimal graphs which differ by a vertical translation 
or a vertical reflection will induce the same area function, so all of these solutions 
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should be regarded as equivalent. In §1, I show that, aside from the obvious case 
of two planes of the same slope, there are essentially three other cases of pairs of 
inequivalent minimal graphs which induce the same area form. All of these turn out 
to be known minimal surfaces, with analytic continuations to complete embedded 
minimal surfaces in space, though the last family appears to have been considered 
only recently [Ka]. 

The second problem is like the first, but for the Lagrangian 

V(2OT)2-(p2+92_i)2' 

which arose in Professor Calabi's treatment of extremal isosystolic metrics on sur- 
faces [Ca2]. Here, the answer is not as complete. I have only been able to show 
that, up to some simple equivalences, there is at most a three parameter family 
of pairs of inequivalent solutions of the Eulcr-Lagrange equation which induce the 
same Lagrangian. Though the analysis could, in principle, be completed without 
the introduction of new ideas, the calculation which must then be carried out is 
beyond the range of any computer to which I have access. The motivation for this 
analysis is explained in the introduction to Section 2. 

The third and final problem to be considered here is of a different nature, though 
it is in the same general class. Given a (simply-connected piece of a) Riemannian 
n-manifold S" of constant sectional curvature K, what are the possible constant 
values that can be assumed by the energy density of a harmonic mapping 

/ : S ^ S'^? I show that, if K is positive, and is constant, then ||(if|P = 

m{m+n—l)K for some integer m > and that, moreover, such an / is, in a suitable 
sense, the restriction to S of harmonic polynomials of degree m on Thus, in 

this setting local harmonic mappings of constant energy density extend to global 
harmonic mappings. This result generalizes results of Wallach [W] and myself [Br]. 
Moreover, the non-uniqueness of / in this case is reduced to a representation- 
theoretic problem analogous to that considered in [DW] and [DZ] in the case of 
minimal isometric immersions. 

1. Minimal graphs with equal areas 

I want to begin by considering a simple, but nonetheless interesting example, 
namely the Lagrangian for the Euclidean area of a graph {x, y) ^ {x, y, u{x, y)), 

L{x,y,z,p,q) = y/l+p^ + q'^. 
The Euler-Lagrange equation for the associated area functional is 

(1) — I I _^ A ( •"y 1-0 

dx y^jYT^^j dy y^j^^-T^j 

Given an open, connected domain T) in the xy-plane, and a (sufficiently differ- 
entiable) function F on 2? which satisfies F >1, the problem is to determine how 
many functions ui5iT> satisfy the minimal surface equation as well as the first order 
PDE 



(2) 



l + ul + ul = F\ 
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If u is one solution of (1-2), then all of the fimctions of the form ±u + c (where c is 
a constant) are also solutions of (1-2), but I will regard these as equivalent to the 
original solution. The problem is to find those functions F for which there exist at 
least two inequivalent solutions, i.e., two minimal surfaces graphed over the same 
domain V in the plane which induce the same area function F onT) without being 
congruent by vertical translation or reflection. 

Before stating the result, let me make a few comments about some evident 
symmetries of the problem. If u is a solution of the minimal surface equation on 
a domain V, then for any constant A > 0, the function ux{x,y) = Xu{x/X,y/X), 
defined on the scaled domain XV is also a solution of the minimal surface equation 
(its graph is the A-dilation in space of the graph of u). Moreover if L{u) = F, then 
L{u\) = F\ where Fx{x,y) = F{x/X, y/X). Thus, if F allows two inequivalent 
solutions, then so will Fx- I will say that the functions Fx are equivalent up to 
scale. A similar remark applies to translations or rotations of the domain V in 
the xy-plane. It obviously suflaces to solve the problem up to these sorts of trivial 
equivalences. 

Theorem 1. Suppose that a function F on a connected plane domain T) is the 
area function of two inequivalent minimal graphs over V. Then, after possibly 
translating, scaling, and rotating in the xy-plane, F is the restriction of one of the 
following functions to a subdomain: 

(1) F = c where c > 1 is a constant andV = M?. There is 1 -parameter family of 
corresponding inequivalent minimal graphs, namely, the circle of planes 
defined by the equations 



(2) F = cotha; and V = {{x,y) \ x > 0}. There is a 1 -parameter family of 

corresponding inequivalent minimal graphs, namely the graphs each of 
which is the part lying over V of Scherk's fifth {minimal) surface 



(3) The function Ffj, on the domain Vcj, = {{x, y)\x'^ + y'^ > sin^ 0} defined by 



There are exactly two inequivalent (multi-valued) minim,al graphs associated 
to F^. They are the upper and lower sheets of the "heli-catenoid" described 
below. 

(4) The function Fa,c defined on T>a.c = {( r. y) \ a cosh.x + c cosy > 1} by 



where the constants a and c satisfy \a—c\ < 1 < a+c. There are exactly 
two inequivalent minimal graphs overV^^c associated to F^^c- Up to vertical 
translation and reflection, they are subsets of a single complete, doubly- 
periodic minimal surface in M?, as described below. 



z = 




sinhx sinhz = cos(j/ + ip). 





The remainder of the section is devoted to the proof of this theorem. 
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1.1. Reduction to an equation. First, a couple of preliminary remarks. If 
a continuously differentiable function u satisfies the minimal surface equation (1) 
on a connected domain V c M."^, then u is known to be real-analytic. If u also 
satisfies (2), then it follows that F must be real analytic as well. Thus, it may 
assumed from the start that F is real-analytic. 

Moreover, I may also suppose that the locus F = 1 consists of isolated points, 
since, otherwise, either F = 1 (in which case it is clear that the only solutions 
to (2) are u = c) or else the set F = 1 would contain an analytic arc a CV. The 
graph of u would then be tangent to a (horizontal) plane along a. However, two 
connected minimal surfaces which are tangent along a curve are identical, so the 
graph of u would be a horizontal plane, implying that F = 1, which, as has been 
seen, does not furnish a solution to the problem. Thus, by excising isolated points 
where F = 1,1 may assume from now on that F > 1 at every point of V. 

Set F = coth/u where is a positive (analytic) function.^ Then, for any solu- 
tion u of (1-2), there must exist a function 9 onV (well-defined up to an additive 
multiple of 27r) so that 

cos 6 , a'mO 

(3) Ux = -T—, — and Uy = — — . 

smh /X smh /z 

Equation (1) becomes 

d / cos6i \ 9 / sin6' \ _ ^ 
dx \ cosh n J dy \ cosh ^ J ' 

and the identity {ux)y = {uy)x becomes 

sin^ \ 9 / cos 6* 



= 0. 



dx \sinh/i/ dy \sinh/i^ 
Solving this latter pair of equations for the partials of yields the formulae 

sinh2/z 9x = sin(20) fix — (cosh2/z + cos(2^)) Hy, 
sinh2/x dy = — sin(26') fiy + (cosh2/i — cos(20)) /ix- 

Differentiating the equations (4) and equating mixed partials of 9 yields the 

relation 

(5) {fixx - ^yy) cos(20) + 2fj.xy sin(2^) = {fixx + fJ-yy) cosh(2^) 

which evidently would determine 9 modulo tt up to a choice of at most two solutions 
unless the coefficients of cos(2(?) and sm{29) in (5) were to vanish identically. This 
observation motivates breaking the analysis into two cases. 

1.2. Degenerate cases. Suppose that n satisfies the condition that (5) is an 
identity as an equation for u. This is equivalent to requiring that fj, satisfy 

l^xx l^xy l^yy ^* 

In other words, is a linear function of x and y: 

Ijl{x, y) = ax + by + c 
for some constants a, b, and c (not all zero). 



^Just why the equations that follow are drastically simpler in terms of /j, than they would be 
if expressed directly in terms of F is a mystery to me. 



6 



ROBERT L. BRYANT 



1.2.1. Planes. If a = & = 0, then ^ ~ c\s a, constant, so the equations (4) imply 
that ^ is a constant and hence, by (3), that u is a hnear function of x and y. Of 
course, this is the first case of the theorem. 

1.2.2. Scherk's Surfaces. On the other hand, if at least one of o or 6 is non-zero, 
then by rotating, translating, and scaling if necessary, I can reduce to the case 
H = X. Since /x > 0, it follows that I? is a subdomain in the right half plane a; > 0, 
so I may as well let T) be this entire half-plane. The equations (4) can now be 
integrated, with the result that there is a constant V so that 

6{x, y) = — tan~"^ (tanh(x) tan(j/ + tp)) . 

Now integrating (3), the corresponding 1-parameter family of inequivalent solu- 
tions u is 

u(x,y) = sinh""^ / cos(y + ij. ) \ ^ 
\ smhx ) 

Thus, the graph of u is one-half of one of the classical Scherk surfaces, defined in 
xyz-space by the equation 



sinhx sinh2 = cos(y -|- V')- 



Note that the corresponding area form is dA = coth xdx/\dy and is invariant under 
translation in the ^/-direction. This constitutes the second case of the theorem. 

1.3. Higher compatibility conditions. Now, I want to turn to the case where 
equation (5) is not an identity. In order for there to be any solutions 0, the quantity 

(6) P = {fixx - IJ-yy)'^ + (2Mxj/)^ - C0sh^(2/i) {fixx + IJ-yy)"^ 

must be non-negative. Set A = {ii^x — l^yyY + C^l^xy)^- If A = 0, then P > 
implies that iixx = f^xy = fJ-yy, which has already been studied. Thus, from now on 
I assume that A > 0. The two solutions of (5) are found to be 

/ cos 20 \ / Hxx - l^yy \ (fex + Mw) cosh _^ / -2/i^y \ VP 

\sm2e J \ 2lixy ) A yiJ-xx - IJ-yy ) A 

From now on, I am going to assume that P > 0, since otherwise, there is at most one 
solution modulo integral multiples of tt and, consequently, at most one minimal 
surface (up to equivalence) whose area function is F . Substituting each of these 
two solutions into the equations (4) yields four equations of third order for the 
function /i. These equations arc fairly ugly and I won't write them out. What is 
remarkable is that making a change of dependent variable by setting C = cosh2/z 
results in equations for C which are quite simple: 

Cxxx ~ iS-^x ^xx ^x Cyy ~t~ ^y ^xy) 1^ 
^xxy = (px ^xy)/C 
Cxyy = {Cy Cxy)/C 

Cyyy — Cyy Cy Cxx ~t~ Cx Cxy^ j C . 
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Now, in general, one doesn't expect an overdetermined system such as (8) to have 
any solutions at all, but this system turns out to be involutive (as is easily checked) 
and hence at any point {xq, yo), one can freely specify the values of 0), Cx, Cy, 
Cxx , Cxy , and Cyy and get a unique solution to the above system of equations. It 
follows that there is a six parameter family of solutions. Moreover, these equations 
are seen to have the following first integrals: For any solution of (8), there exist 
constants ai, 02, 03 so that 

C'kj/ = di C 

(9) (^xx — Cyy = a2 C 

C [Cxx + Cyy) — Cl — Cy = as 

Moreover, for any choice of the constants a,, any solution of (9) also solves (8). 

To write out the general solution of (9) , I will take advantage of the invariance of 
these equations under the group generated by translations, rotations, and scaling in 
the a;?/- variables. Rotating the independent variables x and y by an angle (p rotates 
the pair (ai, 02) by an angle of 2(p. Thus, I can assume that ai = and 02 > 0. 

1.3.1. The Heli-Catenoids. If (01,02) = (0,0), then the solutions of (9) are of 
the form 

C{x, y) = a{x^ + y'^) + bx + cy + d 

for some constants a, b, c, and d. It is easy to see that if a = 0, then P = 0, a 
case already discarded. Thus, a 7^ and so, by translation in the .xy-plane, I can 
assume 6 = c = 0. By direct calculation, the condition P > is now seen to imply 
that d^ < 1. Since C = cosh2/Lt > 1, it then follows that a must be positive. By 
scaling in the xy-plane, I can therefore assume that a = 2 and take d = cos2(^. The 
corresponding formula for F is then found to be 



(10) F^{x,y) 




y"^ + cos^ 



x^ +y^ — sin 



where < < 7r/2 and of course, V is defined by the inequality x^ +y^ > sin^ (p. 

A moment's thought leads to a guess as to what the corresponding minimal 
surfaces are: They are the surfaces in the isometric family of associated minimal 
surfaces containing the helicoid (</> = 0) and the catenoid (0 = 7r/2), i.e., the so- 
called "hcli-catcnoids" . The two inequivalent solutions u are got by simply following 
the upper and lower halves of the catenoid through the deformation in the associated 
family. A direct calculation, which I won't do here, verifies this guess. 

1.3.2. Doubly Periodic Solutions. If ai = but 02 > 0, then by scahng in the 
xy-plane, I may further assume that 02 — 1- The equation Cxy — then implies 
that C{x,y) — a{x) + (3{y) for some functions a and (3 of single variables. The 
equation Cxx — Cyy = C becomes 

a"ix)-aix)=f3"{y)+m- 

It follows that each side of this equation must represent a constant function. By 
subtracting this constant from /3 and adding it to a, I can arrange that a"{x) = a{x) 
and l3"{y) = —(3{y). The third equation of (9) is now an identity. 
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The conditions that P > and that C > 1 on 2? then imply that appropriately 
translating in x and y will put the solution in the form 

C{x,y) = a cosh a; + c cosy 

for some constants a and c satisfying the inequalities \a—c\ < 1 < a+c. (This 
describes a scmi-infinitc strip in the first quadrant of the ac-planc.) The corre- 
sponding domain Pa.o i-e-) where C > 1, is the xy- plane minus a sequence of 
congruent 'ovals' centered on the points {x, y) = (O, {2k + l)n). (The size of these 
'ovals' depends on the values of a and c.) 
The function F then takes the form 



(11) Fa,e{x,y) 



I a cosh X + c cos y + 1 
a cosh x + c cos y — 1 



To get the function 9, or more precisely, the functions cos 6* and sin6' (which are 
all that are actually needed), one must now solve (5), which has now become 

(c^ + 2ac cosh a; cosy + — l) cos 20 + 2ac sinha; siny sin 20 

= a{a^ — — 1) cosh a; + c (a^ — + 1) cosy. 

The analysis of this equation is simplified by introducing the smooth surface Ea,c 
in defined by the relation 

z'^ = a cosh a; + c cosy — 1. 

Away from its intersection with the plane z = 0, this surface is two smooth copies 
of 'Da,c, and may be pictured roughly as two copies of the a;y-plane joined by a 
countable number of "throats" evenly spaced along the y-axis. One computes that 
cos 20 and sin 20 can be smoothly and globally defined as functions on Sa,c so as 
to satisfy the above equation by means of the formulae 

cos20\ 1 /A -B\ f E 



sin 20 y A^ + B'^\B A J \Qz ^ 

where A, B, E and Q are smooth functions on the a;y-plane defined by the formulae 

A = + 2ac cosh x cos y + — 1 
B = 2ac sinh x sin y 

E = a{a^ — — 1) cosh a; + c(a^ — + 1) cosy 

Q = ((1 - (a - c)2)((a + c)2 - l)(a coshx + c cosy + 1))^ 

Note that A"^ + B'^ and Q are strictly positive on Sa,c, so these formulae do, indeed 
define cos 20 and sin 20 as smooth functions on T,a,c- By construction, these two 
fmictions arc periodic in the y-dircction, of period 27r, and have the same values at 
the points (a;, y, z) and (—a;, — y, z). 

Now, I claim that it is possible to define cos0 and sin0 as smooth functions 
on J^a,c- To see this, it sufiices to show that, as a multi-valued function on J^a,c, 
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the function changes by an integral multiple of 27r as one traverses any closed 
loop on T,a,c- Because of the j/-periodicity, it suffices to check that this happens 
for two specific loops, say a and 7, where a is any one of the loops on the surface 
cut out by the plane a; = and 7 is any one of the loops on the surface cut out by 
the plane z = Q. For the first loop cr, it is not hard to see that 9 does not change 
value at all when one goes around it. The second loop is somewhat more difficult 
to understand directly, so 1 will replace it by a homotopic loop 7^. Fix a large 
number R ^ Q and consider the curve 7^ on the upper sheet of lla,c which lies 
over the boundary of the rectangle |a;| < i? and < y < 2-k. The corners Pi, P2, 
P3, and P4 lie above the points (— P, 0), (P, 0), (P, 2-k), and (— P, 2it), respectively. 
Now, dividing the original equation above by and assuming x = R (where P is 
large), one sees that 



cos{2e{R,y,z)-y) 



1 



2c 

while dividing the original equation by and setting x = — P, one sees that 

.2 ^2 1 

cos{26{-R,y,z) + y) « 



2c 

where the approximations are uniformly close in P. Since — — 1| < |2c|, it 
follows that when a; ^ the function 9{x,y,z) is approximately ^y + 9q where 
is a constant and when a; ^ the good approximation is —^y + 9o- Using this 
and the evenness of cos2^? and sin 20, one sees that, as one goes counterclockwise 
around the path ^r, the multivalued function 9{x,y,z) increases by 27r. This, plus 
the 27r-periodicity in y of the functions A, B, E, and Q, shows that cos^ and sin^ 
can be smoothly defined on the entire surface Sa,c- Moreover, they are periodic 
in y of period An and satisfy 

cos 9{x, y + 27r, z) = — cos 9{x, y, z) 
sin 9{x, y + 2tt, z) = — sm9{x, y, z) 

Now, the desired function u on 'Da,c satisfies 

cos 9dx + sin 9 dy 



(12) du = 



i(a cosh a; + c cosy — 1) 



where the 1-form on the right hand side of (12), say (, can be shown to be globally 
and smoothly defined on Ea,c- Note that when |a;| is large the 1-form C, is quite 
small, so that the desired function u is nearly constant for \x\ large, i.e., its graph 
is asymptotic to a plane. 
Set 



It is not hard to show that Xa,c is non-zero (which sign it has depends on which 
choice of (cos 9, sin 9) was made). Let Lo,c be the lattice in generated by (0, 47r, 0) 
and (0, 0, Xa,c)- Then the image of T,a,c in ^ /La,c is easily seen to be a torus Ta,c 
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with four punctures. The 1-forms dx, dy, and C arc well-defined on Ta^c and the 
period vectors of the vector-valued 1-form ^ on Ta,c defined by 




lie in La,c by construction. Thus there is a smooth map X : Ta^c ^ K'^/^a.c 
so that dX = ^. By its very construction, the image of X is a minimal surface 
in M^/Lq^c endowed with its quotient flat structure as an abelian group. More- 
over the information derived so far shows that this minimal surface is complete 
and of finite total curvature. In fact, endowing Ta,c with the conformal structure 
and orientation induced on it by its immersion into M? /La,c, one sees that it is 
biholomorphic to a Ricmann surface of genus 1 punctured at four points. 

This analysis can be carried further to show that the Weierstrass representation 
of this minimal surface is such that the completion of Ta^c to a curve of genus 1 can 
be identified with C/A, where A C C is a rank 2 lattice generated by a real and 
a purely imaginary period. It is not hard to see that C, = ^e{2dw) where dw is a 
holomorphic diflferential on C/A. Letting gi: C/A — » 5^ ~ CP^ be the meromorphic 
function which represents the Gauss map of the minimal surface, one sees that 

dx\ I {g- 9~^) dw \ 

dy = Re i (g + g-'^) dw . 
C/ \ 2dw J 

Note that g has two simple poles at two of the half-lattice points and two simple 
zeros at the other two half-lattice points. Thus, these minimal surfaces are of the 
type studied by Karcher [Ka]. 

To recover the two minimal graphs, let Sa,c C denote the inverse image of 
the minimal surface X(Ta,c) C R^/La.c- After a suitable vertical translation, the 
part of Sa.c which lies between the two planes z = and z = Ao,c will double 
cover I?a,c with a fold along the 'oval' boundaries. These two halves will then be 
the two inequivalent minimal graphs which induce the same area function. 

This completes the proof of Theorem 1. 

2. Special extremals for isosystolic metrics 

2.1. The isosystolic ratio. I now want to describe a second problem of this 
sort which is motivated by Professor Calabi's work on extremal isosystolic metrics 
on surfaces. For any Riemannian metric g on an n-manifold M , Gromov [Gr] has 
defined the systole Sys{M,g), to be the length of the shortest closed curve 7 in M 
which is not null-homotopic. A manifold M is said to be essential if there is a 
constant c(M) > so that, for all metrics g on M, 

Sys(M,ff)<c(M) (Vol(M,5))'/". 

The smallest possible value of c{M) is called the minim,um, isosystolic ratio for M. 

While it is true that every compact surface other than the sphere is essential, 
there are only three surfaces for which the extremal isosystolic ratio is known ex- 
actly, the real projective plane, the torus, and the Klein bottle. Professor Cal- 
abi [Ca2] has set himself the task of computing or closely estimating the isosystolic 
ratio for other surfaces and, along the way, has turned up quite an interesting 
problem in the calculus of variations. 
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2.2. Professor Calabi's work. Let mc explain how this calcuhis of variations 
problem arises. In searching for an extremal isosystolic metric, Calabi has narrowed 
the search to a class of (singular) generalized Riemannian metrics on a given sur- 
face M which have the property that, aside from a small singular set, the surface 
can be covered by "tight bands" of shortest non- trivial geodesies, i.e., geodesies 
which are systole-long. Outside of this singular set, the surface can be broken up 
into what Calabi calls "fc- regular domains" Uk through which k such bands pass. 

Calabi first shows that, in any 2-regular domain U2 for an extremal isosystolic 
metric g, the two bands must meet at right angles and that the metric g must be 
flat there. He does this by showing that, otherwise, he can deform the metric in 
a neighborhood of such a point so as not to decrease the lengths of the shortest 
geodesies while at the same time decreasing the area, thus violating the assumption 
that g was extremal. 

Calabi then turns his attention to the problem of writing down the possible 
metrics in a 3-regular domain C/3. He shows that each of the three bands Bi n U3 
is calibrated by a closed 1-form of unit length w*, i.e., a closed 1-form of length 
one with respect to g which restricts to each of the geodesies in the band to be 
the differential of arc-length. Of course, these three 1-forms can be written locally 
in the form = dx'' . We can take the first two as coordinates, say = x and 
x'^ = y, and the consider the third as a function z = F(x,y) of the first two. The 
condition that the metric g make each of dx, dy, and dz be of unit length is easily 
seen to be the condition that 

\adx + l3dy\'^^ = + + 2f{x,y) a/3 

where / is given by the expression 



f{x,y) 

and the corresponding area form is 



1 _ _ ^2 

X y 

2FxFy 



dAg = — ^ ^ dx A dy. 

{2zxZyf-{zl+zl-lf 

It follows that, if the metric g is to be extremal isosystolic, then the function z = 
F{x, y) must be a solution of the Euler-Lagrange equation of the functional 

•^('^) — f — dxAdy, 
l{2zxZyf-{zl + zl-lf 

at least in the region where the three bands are tight. 

It is not difficult to verify that the Euler-Lagrange equation for this functional 
is elliptic at solutions which satisfy the inequalities z1^-Zy > 1 and \z'^—Zy\ < 1. 
Now, because the Lagrangian 

2pq 



L{x,y,z,p,q) = 



V(2pg)2-(p2+52_i)2 
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does not involve x, y, or z explicitly, the Lcgcndre transform of its Euler-Lagrange 
equation becomes an equation for a function of p and q which, on the region £ in 
pq-space defined by the inequalities p'^+q'^ > 1 and < 1, is a second-order, 

linear, elliptic pde of the form A<I> = where A is the Laplacian of an appropriate 
metric on £. Of course, this equation can be reduced to Laplace's equation once 
one finds a complex coordinate ( : 8 ^ C which uniformizes the induced conformal 
structure. (This is the analog of introducing conformal coordinates on a minimal 
surface in space and deriving the Weierstrass representation for minimal surfaces: 
The equation is reduced to a linear one by a change of coordinates.) However, 
such a coordinate has proved rather difficult to find. In fact, no solutions of the 
corresponding Euler-Lagrange equation other than those which lead to a flat met- 
ric g are known, and these trivial solutions are useless for finding a uniformizing 
coordinate on £. 

2.3. The prescribed Lagrangian problem. To make further progress, at least 
one non-trivial solution to the Euler-Lagrange equation needs to be found. The 
usual approach for finding non-trivial solutions would be to look for special solu- 
tions, particularly solutions which are invariant under some continuous group of 
symmetries of the equation. Unfortunately, since the Lagrangian has relatively few 
symmetries, one cannot expect solutions with continuous symmetries and, in fact, 
no non-trivial solutions invariant under a continuous symmetry of the equation are 
known. 

However, returning to the original problem, Calabi realized that understanding 

the extremal isosystolic metrics on 4-rcgular domains involves finding two solutions 
z and w to the Euler-Lagrange equation which induce the same metric g, i.e., which 
satisfy 

2ZxZy _ _|_ 2Wx Wy 

^{2z,zyr -{zl + zl-iy ~ ^{2w,WyY -{wi+wi-iy' 

More generally, in a fc- regular region (where > 3), one expects that for the area 
form dAg = F{x, y) dxAdy there should be fc— 2 distinct solutions z^ to the Euler- 
Lagrange equation which also satisfy 



^(2 44)2 -((4)2 + (4)2 _ 

Thus, this is a prescribed density problem of the sort under discussion. Again, since 
the Lagrangian does not contain z explicitly and because it is an even function of 
{p,q), it follows that, for any solution z to the problem, every function of the 
form ±z -|- c is also a solution. However, such solutions are to be regarded as 
equivalent because they give rise to the same local foliation by geodesies. I can now 
state the main result: 

Theorem 2. For any non-constant function F on a domain V in the plane, there 

are at most two inequivalent solutions to the Euler-Lagrange equation for L which 
satisfy L = ±F. Moreover, up to translation and dilation, the space of those F for 
which there are two inequivalent solutions to this Euler-Lagrange satisfying L = ±F 
is of dimension at most 3. 

This theorem has the following corollary: 



ON EXTREMALS WITH PRESCRIBED LAGRANGIAN DENSITIES 



13 



Corollary. For any extremal isosystolic metric g, the Gauss curvature vanishes 
in all k-regular domains with k > 5. Moreover, there is a family T of metrics, 
depending on at most 3 parameters, so that in a i-regular domain an extremal 
isosystolic metric g is locally isometric to one of the metrics in this family. 

I conjecture that, in fact, the space of F for which there are two incquivalcnt 
solutions consists only of the constants. Of course, if this conjecture is true, then 
the corollary about local flatness can be sharpened to include the 4-regular domains 
as well. This would leave only the 3-regular domains as unknown open subsets of 
the extremal isosystolic metrics. 

The proof of Theorem 2 rests on an analysis of the compatibility conditions for 
the unknown function F with L given as above. These calculations are quite long 
and were carried out on a computer, using the symbolic manipulation program 
MAPLE. For this reason, at a certain point, I am only going to indicate how the 
calculations are to be carried out. 

2.4. Reduction to an equation. Here is the setting of the problem in the elliptic 
case: For a given function F > 1 on a domain V in the xy-plane, find all of the 
inequivalent solutions z on V which satisfy the Euler-Lagrange equations for the 
functional 

C{z) = / — ^ dxhdy 

J-D ^{2zxZyY-{zl + zl-lf 

and which also satisfy the first order equation 

2 Zx Zy 



^{2zxZyf-{zl + zl-iy 



= F. 



Setting F = \/ sin(2(/)) where (/> is a function on V with values in the range 
Q < 4> < 7r/4, the above first order equation can be put in the form 

{2cos2cl,ZxZyf = {zl + zl-lf. 

Considering the elliptic range and reflecting in x and/or y if necessary, I can assume 
that Zx and Zy are positive and satisfy Zx+Zy > 1 and Iz^—Zy] < 1. Then taking 
the positive square root of both sides yields 

zl — 2 cos 2(j)Zx Zy + Zy = 1. 

Thus, (zx,Zy) is determined to lie on an ellipse which passes through the points 
(1,0) and (0,1). From this and the ellipticity assumption, it is easy to see that 
there must be a function 6 satisfying \6\ < ^tt — so that 

(1) Zx= ''''[^:/\ and zy='°'^^ + '^^ 



sin 20 sin 20 

Now, the Euler-Lagrange equation is equivalent to the condition that the 1-form ^ 
be closed where 

^. ^ (4 - 4 - 2zl + l)zx dx - (4 -zj- 2zl + l)zy dy 

((2^.^,)2-(^2 + ^2_l)2)f 
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Substituting the formulae for and Zy into the equations dnl) = d{zx dx + Zy dy) = 
and expanding allows one to solve for the partials of the unknown function 9 in 
terms of the known function (f) and its partials. The actual formulae are rather 
complicated, but for what I am going to describe, it suffices to know that the result 
is an equation of the form 

(2) d{2e) = cos 261 wi + sin 26 uj2 + ^3 

where Wi, W2, and W3 are certain 1-forms which are linear combinations of dx and 
dy with coefficients given by some (rather ugly) rational expressions in sin (j), cos <j), 
(^x) and (j)y (and linear in the latter two variables). 
Differentiating (2) yields 



where Ai. A2, and A3 arc some (also rather ugly) rational expressions in sint/), 
COS0, <^i;, 0a;x, <^a;j/, and (and linear in the last three variables). It follows 
that any solution 6 must satisfy the relation 



(This is the analog of equation (5) in §1.1.) 

Just as in the case of minimal surfaces, there are now two cases to consider: The 
first is that (p might satisfy the second order quasi-linear PDE system 



Just as in the case of minimal surfaces, these three equations can be solved for ^^xi 
(j)xy, and (j)yy, so that this system can be written in the form 

(l>xx = Fl{<j), (px, (i>y), (pxy = F^icj), (px, <Py), <f>yy = Fz{(j), (f)x, (py). 

However, unlike the analogous system for minimal surfaces treated in §1.2, this 
system is not involutivc. In fact, an analysis of these equations using MAPLE shows 
that the only solutions of this system are given by setting <j) equal to a constant. Of 
course, in this case, 6 must also be constant and hence z must be a linear function 
of X and y, yielding the known case where there arc inequivalent solutions. 

The second case is that not all of the Ai vanish, so that (3) becomes a non-trivial 
equation for 6. In fact, clearly there can be at most two solutions 6 of (3) which lie 
in the elliptic range \9\ < ■^n — (j). The condition that there actually be two distinct 
solutions is a set of second order differential inequalities on that it is not helpful 
to write out here. 

Supposing that there arc two distinct solutions, say 9^, they can be written out 
explicitly in terms of </> and its first and second derivatives in a formula analogous 
to equation (7) of §1.3. These expressions are rather large and I won't write them 
out here. It remains to check whether these two solutions of (3) actually satisfy (2), 
i.e., that the following equations hold: 



= d{d{29)) = {Ai cos 29 + A2 sin 26* + A3) da; a dy 



(3) 



Ai cos 29 + A2 sin 29 + A3 = 0. 



Ai=A2 = A3 = 0. 



(4) 



d{2e^) = COS 2(9=^ wi + sin 261=^ 



W2 + 1^3- 
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Substituting the explicit formulae for 9^, the equations (4) expand to a system of 
four third-order equations for (j). In fact, these equations can be solved in terms of 
the third derivatives in the form 



Unfortunately, the functions Gi arc gross and I have not been able to check their 
integrability conditions. It could well be that the only solutions are got by setting 
equal to a constant. In any case, only for a solution (j) of this system (satisfying the 
appropriate second order differential inequalities) do there exist two inequivalent so- 
lutions z to the original problem with prescribed area form F dxAdy = esc 2<p dxhdy. 

What can be said with confidence is that there is at most a six parameter family 
of solutions (f) of (5) (which would be the case if and only if this system were in 
involution, as was the analogous system (8) of §1.3). These solutions reduce to at 
most three after taking into account the fact that the whole problem is invariant 
under translation in x and y and simultaneous dilation in x, y, and z. (This implies 
that for any solution <p(x,y) of (5), all of the functions <p{cx + a,cy + h) are also 
solutions of (5).) This observation completes the proof of Theorem 2. 

Of course, it is possible that a clever change of variables will simplify the equa- 
tions (5) to something manageable, but I haven't been able to find one. My sus- 
picion is that the only solutions are the constants, but it looks as though this will 
be hard to prove. In theory, checking the integrability conditions of (5) is routine, 
but, using MAPLE, it takes more memory than any machine to which I have access. 

3. Harmonic maps of constant energy density between spheres 

3.1. Minimal surfaces and harmonic maps. In this last part, I want to con- 
sider another case where the notion of solutions of Euler-Lagrange equations with 
prescribed Lagrangian makes an interesting appearance. Many years ago. Profes- 
sor Calabi [Caljclassified the minimal 2-spheres of constant Gaussian curvature in 
the standard Euclidean spheres, showing that they were simply the already known 
Boruvka spheres. This work was carried further by Wallach [W] , who showed that 
any minimal submanifold M*^ C 5" of constant positive sectional curvature is a 
(piece of a) global isometric immersion /: {S^,r) 5*" given by harmonic polyno- 
mials of the appropriate degree. 

Now, the condition of being a minimal submanifold of constant positive sectional 
curvature is very strong. Recall that when one endows an n-dimensional minimal 
submanifold with the induced metric, the inclusion mapping becomes a harmonic 
map which is also an isometric embedding. In particular, its energy density is a 
constant n times the volume form. Thus, minimal submanifolds furnish examples 
of harmonic maps with constant energy density. Thus, it makes sense to generalize 
the study of minimal isometric embeddings to that of harmonic maps of constant 
energy density. Both of these subjects have received attention in the literature 
from a nmiiber of authors [Br], [DW], [DZ], [Ke], [L], [PU]. [T]. Note that, since 
harmonic maps are just solutions to the Euler-Lagrange equations obtained by 
taking the energy density as the Lagrangian, this sort of problem fits exactly into 
the sort being discussed. 



(5) 




4^X1 ^xxi ^xyi 4'yv)i 
4^X1 ^yi 4'xxi 4'xyi ^yy^i 
G^icj), 4'x,4'y: 4>xx, 4>xy, 4'yy), 
XI 'Pyi 4'xxi 4'xyi 4yy)- 
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In [Br] it is proved that if E is a surface endowed with a metric of constant 
Gauss curvature K and /: T? ^ is a harmonic map whose energy density \df\^ 
is constant, then K >Q. Moreover, \i K > Q then /(E) is a piece of a Boruvka 
sphere.^ Recently, I have extended this result to all dimensions. 

Theorem 3. If {Y7^,K) is a (piece of a) space form of constant sectional curva- 
ture K > and f:T,—^ is a harmonic mapping with constant energy density, 
then f extends to a global harmonic map from 5" to which is given by harmonic 
polynomials on 5" of the appropriate degree. 

The proof, which is based on an interesting identity, will be outlined in the 
following sections. Because the notation used below has certain failings when n = 2, 
I am going to assume from now on that n is at least 3. In any case, the theorem is 
already proved in [Br] in the case n = 2. 

3.2. A Notation for Representations of the Orthogonal Group. The proof 
of Theorem 3 will involve working with symmetric tensors of rather high order on 
a space form E" of dimension n. For this reason, it is important to establish an 

effective notation for the irreducible components of the representations of SO(n) 
on the symmetric powers of M" and various bilinear pairings between them. 

Let x^,. . . ,a;" denote the standard linear coordinate functions on M". Let V 
denote the ring of polynomials in the variables and let Tid C V denote the 
vector space of harmonic polynomials of degree d. It is known that Tid is an 
irreducible representation of SO(n). Moreover, the symmetric tensor power S'^{Hi) 
of Hi ~ M", regarded as the space of all homogeneous polynomials of degree d 
on M", obeys the inductive decomposition rule 

where R= {x^Y^ [-{x^Y, thus implying the complete irreducible decomposition 

of this symmetric power as an SO(n)-module. 

For any / G Tid and any G Hi, the product ^/ is a homogeneous polynomial 
of degree d+1 and hence can be decomposed into irreducible components by the 
above rule. In fact, it turns out to be best for normalization purposes to define 

(n + 2d-2) e/ = /Ve + i?^(/-0 

where / V ^ lies in Hd+i and / • ^ turns out to lie in Hd-i- In fact, / • ^ turns out to 
simply be the directional derivative of / in the direction ^ (regarded, via the metric, 

as a vector in M"). This imambiguously defines two pairings, V: Hd x Hi ^ Hd+i 
and •: Hd x Hi Hd-i, which will be used extensively. For example, when n = 3, 
one obtains 

= 2{x^f - (x'^f - (x^)^ and x^ ■ x^ = 1. 

The Lie algebra so(n) of SO(n) is isomorphic to A2(K") ~ A'^{Hi), and the 
formula for the derived action of this algebra on the representation Ha (regardless 
of the value of d) is found to be 

(QA/3)./ = a(/-/3)-/3(/-a). 

^This rules out there being any minimal surfaces in any S" of constant negative curvature, 
which is reminiscent of a very early result of Professor Calabi's to the effect that the Poincare 
disk cannot be isometrically embedded as a complex curve (even locally) in any CP". 
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Moreover, there is an SO(n)-invariant quadratic form on Tid defined by the rule 



if, 9) 



where A is the (geometer's) Laplacian, i.e., the negative of the trace of the Hessian. 
Finally, one can define an SO(n)-equivariant pairing {, } : Hd x Tid+i Ti-i by the 
rule 

{/, 9}-Oi = {f,g-a). 

A rather large number of identities attend these pairings and the inner product. 
For example, for all / S Hd and a,P GHi, 

{fy a) ■ P - {fy 13) ■ a= {n+2d) (a a f3) . f, 

{f-a)yi3-{f-l3)Va= -(n+2d - 4) (a a /3) . /, 
{fVa)-a-{f-a)Va= (n+2d - 2) | |q| 1^ /. 

Moreover, for all / e Hd and g G Hd+i, the formula 

(/Va, g) = {n + 2d-2) {f,g-a) 

shows that, up to certain well-defined scale factors, the operators Va and -a are 
adjoint with respect to these inner products. Listing all of the identities to be used 
below would result in a rather large, incomprehensible table. Since the proofs of 
these and other like identities are straightforward, I will merely point out when 
they are used and leave the verifications to the reader. 

3.3. Structure equations on a space form. Now consider a local orthonormal 

coframing on (E",ds^), say 

ds^ = {uif + ■ ■ ■ + {u)nf 

By the first structure equation of E. Cartan, there exist unique 1-forms ipij = —ifji 
representing the Levi-Civita connection on the domain of the coframing so that 

doJi = —ifij AUJj . 

(Note the use, here and below, of the summation convention, even though all of the 

indices have been lowered.) Since the metric ds"^ has constant sectional curvature K, 
the second structure structure equation of E. Cartan becomes 

difij = —ipik A ipkj + K Ui AUj . 

Now, these equations can be written in vector form as follows. Let lo and ip 
denote the 1-forms with values in Hi and A'^{Hi) ^ so(n), respectively, defined by 
the formulae 

CO = Ui and <P = 5 ^ij ^ ■ 

Then the first and second structure equations are expressible in the form 

du = — ip . u> and dip = — ^[(p,(p\ + ^ ui aoj. 
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3.4. An identity for eigenfunctions of the Laplacian. Now let / be a smooth 
function on the domain of the coframing. There exists a unique Hi-valued func- 
tion /(^^ on this domain so that 

df = /(I) • u. 

(The components of J*^^) arc basically the components of V/ in the orthonormal 
coframing.) Now, it is not difficult to see that there is a formula of the form 

d/W= -^./a)--^A/Va; + /(2).a;, 

where A/ is the (geometer's) Laplacian and /^^^ is a W2-valued function which, 
up to scale, gives the components of the traccless part of the Hessian of / in this 
coframing. This formula could be supplemented by formulae for df'^^^ etc, but these 
rapidly become complicated for a general function. However, for eigenfunctions of 
the Laplace operator, there is a very nice sequence of formulae: 

Lemma. Let f satisfy A/ = A/ for some constant A. Then there exists a sequence 
of functions f = /(°), /^^^, f^^\ with Z*^™) taking values in Hm, so that 
= /(I) . u and 

(1) d/W=-^./W-6™_i(A)/('"-i)Va; + /("+i).a; (m > 1), 

where the seQUCTice 6771, (A) for tti ^ is given by the fovTuuld 

A — ra{n -\- m — 1)K 



(2) bm{X) 



(n + 2m)(n + 2m-2)' 



Moreover, the function f satisfies y("»+i) = Q if and only if it is the restriction 
to T, of a homogeneous harmonic polynomial of degree m. [when S (or its simply 
connected cover) is given its usual interpretation as a hypersphere in a vector space 
of dimension n+1). 

The proof of the formula part of the Lemma is a straightforward induction, 
using the pairing identities listed above and a few others very like them. I will not 
reproduce it here. Of course, the function /^"^^ is, up to a scale factor, the traceless 
part of the fully symmetrized m-th covariant derivative of /. 

The formula (1) can be partly interpreted as a formula relating certain first order 
differential operators 

V±:C°°(S'^(T*)) — > C°°(5'^±i(T*)). 

In this language, suitably interpreted, equation (1) implies that if / € C°°(E) 
satisfies V^V+/ = A/, then 

V-(V+)'"/ = 6„-i(A) {V+r-'f. 

However, the Lemma contains more information than this. Moreover, the defi- 
nitions of the operators V*^ depend on choices of various scale factors which are 
notoriously inconsistent in the reference literature. Thus, I have chosen to avoid 
this way of explaining the proof 
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Seeing the second part of the Lemma, about the identification of the eigenfunc- 
tions which satisfy = 0, require a few remarks. First, of all, notice that, 

by (1), if m is the lowest value for which vanishes, then 6m(A)/^™^ must 

vanish identically, so, by the minimality of m, it follows that bm{^) must vanish, 
i.e., that 

A = m(m + m — 1)K, 

which is exactly the eigenvalue of the restriction to E of a harmonic polynomial 
of degree m on the ambient vector space in which S is embedded in the standard 
embedding. Moreover, setting = and considering the remaining differ- 

ential formulae for /(°\ . . . , /*^™', one sees that this becomes a set of linear total 
differential equations which is closed under exterior derivative. It follows that the 
vector space of solutions to this system is of dimension 

dimWo + dimWi + • • • + dimW„ = ^„ 

m!(n — Ij! 

which, again, is exactly the dimension of the space of harmonic polynomials on 
jjn+i q£ (jQgj-QQ Of course, this information is not conclusive. The remainder of 
the argument comes from the following formula, which is easily derived from the 
differential identity of the lemma and various pairing identities: 

^l^^^im)j{m+i)^, *a;) 

= (m + l)(/(™+i),/('"+i)) -6„(A)(n + m-2)(n + 2m)(/('"),/('"))l *1. 

(Here, * is the Hodge star operator.) Because of the SO(n)-invariance of all of 
the quantities involved, the n-forms on either side of this equation are well-defined 

on S, independent of any choice of local coframe. Thus, if S is compact, then 
integrating both sides of this equation over S yields 

I ^j,(m+l) j(m+l)^ _ bm{X){m + 71- 2)(n + 2m) f ^^(m) jim)-^ 

Ji: (™ + 1) Js 

Prom this, it follows immediately that if A = m(rn + n — 1)K , then any global 
eigenfunction / of the Laplacian with this eigenvalue (in particular, the restrictions 
of the harmonic polynomials of degree m on IR"'+^) must satisfy {f^"^~^^\ = 
and hence = 0. The remainder of the Lemma now follows by dimension 

count. 

It is amusing to note that the above integral formula shows directly that there 
cannot be any non-zero eigenfunctions of the Laplacian on a compact S for any 
value of A not of the form m(m-|-n— l)ii'. The reason is that, if A is not of this form, 
then the fact that K > {) implies that there will be a first value of m for which 6m (A) 
is negative. Of course, this will imply that the integral on the left hand side of this 
equation is not positive, forcing /('"+!) = 0, which, as has been shown, cannot 
happen for any non-zero solution of A/ = A/ unless A = m(m + n— 1)K. 

3.5. Harmonic maps of constant energy density. Now suppose that /: S 

is a harmonic mapping with constant energy density. I will use (, ) to denote 
the inner product on M"+^ which renders the unit sphere.. Thus, / can be 
regarded as an M^+^-valued function which satisfies (/, /) = 1. 
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By the usual computation for harmonic maps to spheres, the condition that / 
be harmonic as a map from E to is 

A/ = IM/||V 

where is the energy-density of the mapping and A is the (geometer's) Lapla- 

cian. Thus, the condition that the energy-density be constant is equivalent to the 
condition that the components of / be eigenfunctions of the Laplacian, all with the 
same constant eigenvalue, say, A. 

Looked at another way, one can say that constructing a (local) harmonic map 
of constant energy density from S to some is equivalent to choosing 7V+1 
(local) solutions of the equation A/o = A/o with the property that the sum of 
their squares is a constant function. Essentially, this can be thought of as A''+l 
second order pde for TV unknown functions, i.e., a (mildly) overdetermined partial 
differential equation, the sort that the theory of exterior differential systems [BCG] 
was designed to treat. 

However, in this case, the problem is more easily treated directly. In the formulae 
below, the notation ((, )) will denote the tensor product of the inner products (, ) 
on Hd and (, ) on ]R^+^. Similarly, given a function g with values in (g) Hd 

and h with values in R^"'"^ (^Hd+i, the expression {{g,h)} will denote the function 
with values in Hi which is got by applying the inner product (, ) and the pairing 
{, } defined above in the obvious way. 

Now, because / satisfies the vector equation A/ = A/, the formulae derived in 
the last section apply directly with little modification to yield, for all m > 0, 



(m + l)((/('"+^\ /('"+^))) - bm{X){n + m-2){n + 2m)((/W, /^^ 



*1. 



Now consider the identity (/, /) = ((/*•"■', /*•"■')) = 1. By differentiation, this implies 
that {(/^°\ Z*-^^)} = 0. Substituted into the above formula, this recovers the already 
established identity ((/^^\ Z'-^-')) = A. However, now differentiating this identity 
yields {(/^^\ /^^^)} = 0. Setting m = 1 in the above relation now gives 

y(2) j(2)| ^ 6i(A)n(n + 4) ^^^(i) ^(i,^^ ^ Abi(A) n(n + 4) 
3 3 

so that, in particular, ((/^^\/^^^)) is also constant. 

Continuing now by induction, one sees that starting with the knowledge that 
(^(^fim)j(m)^^ is a constant function, say, ((/('"),/('"))) = A„(A) > 0, one can 
differentiate, yielding {{f^"^\f^"^~^^^)} = 0. Substituting this into the formula above 
yields 

0(m+l) j(m+l)-j^ ^ bm(A)(TO + n - 2){n + 2m) i^jjm) j(m.)^ 

_ hm{\){m + n-2){n + 2m) ^ 

— 7 — T ^m(.Aj. 

(m + 1) 

Thus, each of the functions ((/*•'"•',/'•'"•')) is a constant ^^(A). These constants 
form a non- negative sequence defined inductively by Aq{\) = and 

bm{X){m + n-2){n + 2m) 
^m+i(A) = ^^^^ A^X). 
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Of course, since K is positive, the seqiience 6m (A) must eventually become neg- 
ative. If A were not of the form m{n+m—l)K , then this sequence would become 
negative without ever becoming zero, which would force the sequence Am{\) to 
become negative at some point. This is manifestly impossible, so there must exist 
an m so that A = m{n+m—l)K. Of course, this imphes 

= ^„+l(A) = 

which in turn implies that = Q. By the Lemma, this implies that / is the 

restriction to E of a harmonic polynomial of degree m, so that / extends to the 
entire n-sphere, as desired. This completes the proof of Theorem 3. 

3.6. Non-uniqueness. For simplicity, I will now restrict to the case K = 1 

and assume that S" = 5*". In [Br], a strong uniqueness theorem was proved: 
If / : 5^ — + is harmonic and of constant energy density, then / is actually 
a minimal isometric immersion up to a constant scale factor, thus implying that 
f{S'^) is a Boruvka sphere. For n greater than 2, and TV and m sufficiently large, 
this sort of uniqueness fails, as can already be seen in [DW] and [DZ]. In fact, the 
non-uniqueness in the present case of harmonic mappings follows from the same 
argument: As has been shown, / : 5" ^ is harmonic of constant energy 
density A = m{m+n—l) if and only if / is the restriction to S'" of a mapping 
F : M"+^ — > whose components are harmonic polynomials of degree m. 

In fact, choosing a sequence F",... of functions in Wm(M"+i), they are the 
components of such an F if and only if 

{Fy + {F^f + ■■■ + {F^f = {{xy + ■■■ + = i?™ 

Let D = dim7im(R"+^) and let {/iq:|1 < a < D} be an orthonormal basis for 
with respect to the inner product defined in §3.2. Given any symmetric 
D-hy-D matrix G = {g"'^), the expression h{G) = g°''^ hahp is a homogeneous 
polynomial of degree 2to on R"^'^. Note that, by equivariance, the identity matrix 
Go = ((5"'^) yields h{Go) = cRJ^ where c > is some constant. Since, for large 
enough m, one has 

dimS\W+^) = ('" + 2'^') « dim52(Hd(M"+i)) « 

V / ((n — 1)!) 

(this strong inequality is where n > 2 is used) it follows that the linear mapping 
G h(G) cannot be injective. The inverse image Q = /i^^^i?"'), which is a non- 
empty afEne subspace since it contains c~^Go, is therefore of positive dimension 
and contains a closed convex set C with non-empty interior consisting of positive 
semi-definite matrices. 

Each matrix G € Q gives rise to a harmonic mapping of constant energy density 
as follows. Write G = (.9"*') = where s is also positive semi-definite matrix (of 
the same rank as G. Then the functions = s^^hp satisfy the above quadratic 
relation and furnish a harmonic map f : ^ 5"""^ where 5"""^ C S^~^ is the 
intersection of the unit sphere in with the r-dimensional range of G. 

Again, since n > 2, the dimension of C grows rather rapidly with to, so that it 
quickly exceeds the dimension of SO(n-l-l), the group of domain symmetries of the 
problem. At that point, non- uniqueness of / up to the obvious symmetries of the 
problem is assured. 
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3.7. Non-positive sectional curvature. In [Br], where the n = 2 case of Theo- 
rem 3 was treated, formulae like (1) were used in the case if = to prove that any 
harmonic mapping of constant energy density from an open set in the Euclidean 
plane to an A'^-sphere was the restriction of a harmonic mapping defined on the 
entire plane and was, in fact, given by trigonometric polynomials. It would be 
interesting to know if this result holds in higher dimensions. 

In the case K < 0, it was shown that there were no harmonic maps of con- 
stant energy density from E to for any A'', a result which is also unknown for 
higher dimensions. This result depended in a crucial way on the growth of the 
sequence Am{^) and other sequences which are closely related. It is still possible 
that the result for higher dimensions will follow from such growth estimates, but I 
do not see how to carry them out. 
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